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Abstract 

The problem of the Darrieus-Landau instabihty at a discontinuous deflagration front in a com- 
pressible flow is solved. Numerous previous attempts to solve this problem suffered from the deficit 
of boundary conditions. Here, the required additional boundary condition is derived rigorously tak- 
ing into account the internal structure of the front. The derived condition implies a constant mass 
flux at the front; it reduces to the classical Darrieus-Landau condition in the limit of an incom- 
pressible flow. It is demonstrated that in general the solution to the problem depends on the type 
of energy source in the flow. In the common case of a strongly localized source, compression effects 
make the Darrieus-Landau instability considerably weaker. Particularly, the instability growth rate 
is reduced for laser ablation in comparison with the classical incompressible case. The instability 
disappears completely in the Chapman-Jouguet regime of ultimately fast deflagration. 
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I. INTRODUCTION 



The Darrieus-Landau (DL) instability is one of the most fundamental and important 
instabilities in hydrodynamics [1^. This instability develops at a deflagration front, i.e., a 
front of energy release propagating subsonically due to thermal conduction. Among the 
many different examples of deflagrations we mention chemical flames l|, [3, 1^, sl, ^, 
. aser ab 

3,y,y, 



ation in inertial confinement fusion and in manipulation of nanostructured surfaces 
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12 



, thermonuclear deflagration in supernovae Q, 0, Q, 



19|, as 



well as waves of phase transition. If we perturb an initially planar deflagration front, the 
perturbations grow because of the DL instability bending the front. At the nonlinear stage, 
the DL instability increases the deflagration velocity; recent reviews on the DL instability in 
flames may be found in, e.g., Refs. jsl,]^. The classical DL theory considers an infinitesimally 
thin deflagration front in an incompressible flow Ij. To the best of our knowledge, the first 



„ to =tudy the Muence of ga. co,np.e.io„ on the DL instabiUty was ,nade in Ref. 
Reference [20| treated compression effects as small corrections to the classical DL theory. 
This is reasonable for chemical flames, since even relatively fast laboratory flames have 
the Mach number much smaller than unity both in the fuel mixture and in the burnt gas. 
However, there is an important example of deflagration with intrinsically strong compression 



effects; this is laser ablation in inertial confined fusion [21|, |22|, |23|. In laser ablation, the 
isothermal Mach number is one at the critical surface, where the laser radiation is partly 
absorbed and partly reflected by the fusion plasma. Due to this, the DL instability in laser 
ablation cannot be studied assuming small compression effects. 

In inertial confined fusion, the DL instability typically appears together with the Rayleigh- 



n 



Taylor (RT) instability (7|,l8|,|9|, 



11 



, 1^ . Quite often the DL instability is overwhelmed 
by the RT instability because of the large acceleration of the plasma targets. Still, taking 
sufficiently large targets, one obtains a small acceleration, which tends to zero for an in- 
finitely large target. In that case the RT instability becomes weak, and one can observe 
the DL instability of ablation flow in a pure form. The DL instability in ablation flow was 
encountered in Ref. Sj within the model of incompressible flow (see Eq. (54) of the paper). 
However, Ref. Q| did not state this finding openly; for the first time it was done much later 



in Ref. 



24j. Since then, there has been much interest to the DL instability in inertial con 



fined fusion both in the theory 
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12 



24 



26 



27( 1 and numerical simulations 28l. 1291. l30| 
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In the theoretical papers 



m 



12 



271] the incompressible approximation was used, which 



may be considered only as a qualitative model for the ablation flow. The influence of com- 
pression effects on the DL instability in laser ablation has been considered only recently 



24 



25j, where, unfortunately, the plasma compression was treated only as second-order 



corrections to the incompressible solution. Moreover, the second-order corrections in Refs. 
24 . 25I concerned only terms proportional to the deflagration thickness; this thickness is in- 



finitesimal in the classical DL theory. Such an approach can therefore not describe properly 
the DL instability in laser ablation. One should expect considerable modifications of the 
DL theory because of the flow compression already for a discontinuous deflagration front. 



26| did not 



The case of strong compression was considered in Ref. j26|. However, Ref. 
resolve properly the deflcit of boundary conditions at the deflagration front and the solution 
was built on a hidden arbitrary assumption (see below). Thus, the fundamental problem of 
the DL instability in fast deflagration and laser ablation in a strongly compressible flow has 
remained without a solution so far. 

The above problems with the DL theory in a compressible flow are not surprising, cf. the 
RT instability in inertial conflned fusion. It is well-known that the model of a discontinuous 
ablation front used to study the RT instability contains a deflcit in the boundary conditions 
0, y, y, 0, 11, Q, [3]: the number of unknown values exceeds the number of conservation 



laws at the front by one. This deflcit was obtained flrst in Ref. [7|; the way to overcome 
the trouble correctly was suggested in Ref. [sl within the incompressible model. In the 
case of incompressible flow the extra condition may be easily guessed: it is the condition of 
a constant flame velocity with respect to the fuel mixture jj]. This is the well-known DL 
condition; it may be also proved rigorously for an isobaric flow, see e.g. Refs. js, 4, 5]. The 
problem of the additional boundary condition becomes difficult again as we take compression 
effects into account. 

In order to overcome the trouble, Ref. 201] assumed an extra condition identical to 



the incompressible DL condition (below we show that this was not the only assumption of 
Ref. {20!]). As a result, Ref. 20] obtained that compression effects make the DL instability 



stronger in the model of a discontinuous front. The papers 
as Ref. 



24 



251 ] used the same assumption 



20| in order to study the DL instability in an ablation flow. Still, unlike Ref. |20| . 



in Refs. 



24 



251 ] no influence of the compression effects for the discontinuous de flag ration 



front was found. Reference 2^ demonstrated that the assumption of Refs. 



24 



25| is 
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incorrect: it contradicts the conservation laws at the deflagration front. In order to resolve 



the problem, Ref. 



261 ] took into account energy conservation at the front, which provided 



one extra condition. Still, in that case one obtains also one extra perturbation mode (the 
entropy mode), which means one extra unknown value. Like before, the number of unknowns 
remains larger than the number of equations by one. The solution proposed in Ref. 26] 



L In order to overcome the deficit and 
261] simply omitted the entropy mode 



reproduces the basic elements of reasoning of Ref. [7| 
to obtain a solution to the spectral problem, Ref. 
without any further explanation. This assurnption may be treated as a hidden additional 
boundary condition in the analysis of Ref. 26[. Thus, we come to the question: what is 
the correct additional boundary condition at a fast deflagration front, taking into account a 
possibly strong compression of the plasma or gas. Of course, this condition should go over 
to the DL condition in the case of incompressible flow. When this question is answered, we 
face the next issue: what is the influence of plasma compression on the DL instability in 
aser ablation and fast deflagration? Does it increase the instability, as suggested in Ref. 

or leave it unchanged like proposed in Refs. [2J, |25|]? Or may be we have the 
third option, and compression effects make the instability weaker. The solution presented 
in Ref. [26] suggested a combination of these tendencies. These questions are especially 
interesting for the Chapman- Jouguet regime of fast deflagration, for which the velocity of 
hot plasma or gas is equal to the local sound speed, isothermal or adiabatic. 

In the present paper we solve the problem of the DL instability at a discontinuous defla- 
gration front in a compressible flow. Taking into account the internal structure of the front, 
we derive the additional boundary condition, which is missing in the model of a discontin- 
uous front. The derived condition implies constant mass flux at the front; it reduces to the 
classical Darrieus-Landau condition in the limit of incompressible flow. We demonstrate 
that the solution to the problem in general depends on the type of energy source. In the 
case of a strongly localized source, compression effects make the DL instability considerably 
weaker. In particular, the DL instability growth rate is reduced for laser ablation in com- 
parison with the classical incompressible case. The instability disappears completely in the 
adiabatic Chapman- Jouguet regime of ultimately fast deflagration. 
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11. BASIC EQUATIONS FOR THE DL INSTABILITY 



We start with basic hydrodynamic equations describing a plasma (gas) flow in the absence 
of gravity 

|^ + V-(pu)=0, (1) 

p^ + (pu- V)u + VP = 0, (2) 
(99 

pT— + pTu • VS - V • (kVT) = (3) 
and the equation of state of an ideal gas 

P = ^^CppT. (4) 
7 

We take the equation of energy transfer in the form describing variations of entropy 

^ = Cvln(^), (5) 

where Cp, Cy are heat capacities at constant volume and pressure, 7 = CpjCy is the 
adiabatic exponent and k is the coefficient of thermal conduction. In general, kinetics of 
energy release in a deflagration flow may involve additional differential equations speciflc for 
a particular type of deflagration. Since it is impossible to consider all particular cases in 
one paper, here we take energy release described by some function f2p(p, T). Most of the 
results of the present work do not depend on ^Ir{p, T); we only demand that energy release 
is strongly localized in a narrow zone inside the deflagration front, which is the usual case. 
This demand concerns not only the function Qpt^ but also the derivatives dQp/dp, dQ^/dT. 
In the case of laser ablation the energy release is typically pres ented by (5-function, which 
implies sufficiently strong localization, e.g. see Refs. [3, 0,1211. We also assume that the 
function allows a planar stationary solution consisting of two uniform flows of cold heavy 
plasma (label "a") and hot light plasma (label "c") separated by a transitional region, which 
is the deflagration front. The labels "a" and "c" originate from the ablation and critical 
surfaces in the laser deflagration. Typical internal structure of the deflagration front is 
illustrated in Fig. [H Strictly speaking, laser ablation requires the flow of hot light plasma 



in the form of a rarefaction wave, not a uniform flow 2l|]. Still, by assuming a uniform 
flow, we may consider all possible values of the Mach number in the light plasma from zero 
to unity, which would be impossible with a rarefaction wave. Besides, the model of two 



uniform flows separated by the transition region is quite common in the studies of the DL 



and RT instabihties in laser ablation, e.g. see Refs. 



24 



25 



26| 



The planar stationary deflagration may be described by the integrals 

pUz = PaUa = PcUc, (6) 
P + pul = Pa + PaUl = Pc + PcUl (7) 

One of the main dimensionless parameters in the problem is the expansion factor 

In the case of ablation flow, the laser frequency determines the critical density and the 
expansion factor. The other important parameter is the Mach number in the light plasma 
(gas) corresponding to the adiabatic sound 

Mai = (9) 

In laser ablation, the isothermal Mach number is equal unity pcU^/Pc = 1 in the light 
plasma, and we have the adiabatic Mach number Ma^ = 1/7. Still, in the present work we 
consider a general case of an arbitrary Mach number between zero and unity. The Mach 
number in the heavy plasma follows from (I7j) as 

9 Mai , ^ 

= 9 + ,(9-l)Mar (^"l 

The internal structure of the deflagration front obeys the stationary equation of energy 
transfer 

PaUaT— - — (K—)=nR. 11 

dz dz \ dz J 

Characteristic width of the deflagration front is determined by thermal conduction in the 
hot plasma with typical definition 

Small perturbations (p to the stationary solution '^{z) have the general form 

(p{x, z,t) = (p{z) eyi\i{at + ikx), (13) 
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where k is the perturbation wave number and a is the instabihty growth rate. In the problem 
of the DL instabihty a is a real positive value. Then linearized system ([I]) - (jH) is 

o-p + ^(pwz + puz) + ikpu,j: = 0, (14) 

apUx + pUz—r— + 'ikP = 0, (15) 
dz 

duz /- ~ ^duz dP , . 

apu^ + pu^— + [pUz + pUz)— + -1- = 0' (16) 

apTS + [pUz + pUz)T— + PaUaT— + PaUaT — 

dz dz dz 

- — (.T) + A:.T=^p+^T, (17) 

P = ^^Cp{pT + pf), (18) 
7 

s-ap-/). (19) 

Since the instability develops at the front, then solution to (fT4|) - (fT9|) should decay at infinity 
{z ±oo) in the uniform flows of heavy and light plasma. All coefficients of Eqs. (1141) - 
( TT^ are constant in the uniform flows, where the perturbations take the form 

(^(x, z, t) = (^exp(cTt + i/cx + /iz), (20) 

with /i > in the heavy plasma z —>■ — oo and /i < in the light plasma z —>■ oo. Then, in 
the uniform flows, the system ([T3j) - ( |T9l) reduces to 

ap + p{pUz + pUz) + ikpUx = 0, (21) 

apUx + pUzPUx + ikP = 0, (22) 

CrpUz + PUzPUz + yU-P = 0, (23) 

cjpT~S + p„f/„T/x5 - p^kT + A;2«:T = + ^T. (24) 

In the case of a strongly localized energy source we have dQ^/dp = 0, dQ^/dT = in the 
uniform flows, and the right-hand side in Eq. (p^ is simply zero. In the present work we are 
interested mainly in energy sources of this type. Still, another situation is also possible. For 
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example, Refs. [3l|, l32] considered a first-order reaction of the Arrhenius type. In that case 
perturbations of the energy source are not zero in the hot gas; they may influence structure 
of acoustic modes and modify the whole solution to the problem. As we show below, such 
modifications take place even in the model of a discontinuous defiagration front. Then, in 
order to solve the problem, one has to specify particular kinetics of energy release, which is 
beyond the scope of the present work. 

III. THE MODEL OF A DISCONTINUOUS DEFLAGRATION FRONT 

A. Solution in the uniform flows 

In the present work we solve the stability problem within the classical DL model of a 
discontinuous defiagration front, kL^ << 1. The traditional scaling for the DL instability 
growth rate is 

a = VUak. (25) 

Our purpose is to find the coefficient F as a function of the expansion factor G and the 
Mach number in the hot plasma Mac- Within the discontinuity model we have to solve 
only Eqs. f|T^ - (fT6!) . The equation of energy transfer Eq. (ITTI) is out of the model, since 
it contains thermal conduction and energy release "hidden" inside the infinitesimally thin 
front. Still, as we will see below, the discontinuity model involves uncertainties, which make 
the solution ambiguous and which may be removed correctly only by taking into account Eq. 
f|T7|) . First of the uncertainties is the relation between pressure and density perturbations 
in the so-called "sound" modes in the uniform fiows. Previous works on the subject Refs. 



24 



25 



261] assumed the adiabatic relation 

P = j-p = clp. (26) 
P 

As we can find from Eq. i^^, this is indeed the case for a strongly localized energy source 
with dQfi/dp = 0, dQ^/dT = in the uniform fiows. The power exponent fi in Eqs. ( !20ll 
- (l2l|) scales for the sound modes as /i cx fc, see below. Within the hmit of discontinuous 
defiagration front, the terms with heat conduction in Eq. fl24p become as small as kL^ << 1, 
and the equation describes drift of entropy perturbations: 

pT{a + u,fi)S = 0. (27) 
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The combination {a + Uzfi) stands for drift with the flow, and it is non-zero for the sound 
modes (see the calculations below). Then Eq. fl271) leads to zero perturbations of entropy in 
the sound modes in agreement with Eq. (^E^. We stress that this is not a general case, but it 
holds only for a strongly localized energy source. In the opposite case of non-zero derivatives 
dilfi/dp, dQfi/dT in f l24l) . the perturbations of energy release become dominating in fl24p 
in the uniform flow of hot plasma. In that case the relation of P and p is determined by 
a particular structure of the energy source. We repeat one more time that in the present 
work we are interested mainly in strongly localized sources with dQ^/dp = 0, dQfi/dT = 
in the uniform flows, which leads to the adiabatic relation (25) in the sound modes. 

With these results in mind, we can obtain an equation for the power exponent /i from 



{a + u,p) [{a + u,pf + cl{e - p'')] = 0. (28) 

Solving (|28|1 . we find possible perturbation modes in the uniform flows describing vorticity 
drift 

Pv = -^< 0, (29) 

and sound 



ht plasma, respectively. Mark that Eq. fl30|) is different from the respective 



^"'^ ~ ~UZ^-Ma^ ^ V t5^(l-Ma2)2 + l-Ma^' 
with Pa > and Pc < for " +" and " — " in ( 130|) corresponding to perturbations in the cold 
heavy and hot li 

results in Ref. 20|, l25||, which were written as expansion of (l30l) in powers of small Mach 
number. Ma << 1. Still, the limit of Ma << 1 does not hold for the ablation flow. The 
structure of the perturbation modes in the uniform plasma follows from Eqs. (!22l) - ( !24l) as: 

Pa = -PaUa [jf- + 1 ) M.a, (31) 
\'J apa J 

k 

for the sound mode in the cold plasma; 

Ps = -pcUc (jf- + 1 ) u,s, (33) 

Uxs = i—Uzs, (34) 
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for the sound mode in the hot plasma; and 

Pv 

UxV 



a _ 



(35) 
(36) 



for the vorticity mode in the hot plasma. Density perturbations are related to pressure 
perturbations by Eq. ( !26l) . Equations ( !33l) - ( l36l) in the hot plasma may be also reduced to 
one condition 



k ^ a 



a 



+ 1 



Pr 



+ 



a 



pJJc Uck 



(37) 



for = U^s + UzV, Uxc = Uxs + UxV, Pc = Ps- 

B. Conditions at the discontinuous deflagration front 

The solution in the uniform flows fl3T]) - fl36|) has to be matched at the perturbed de- 
flagration front Zf = f exp{ikx + at) using the conservation laws of mass and momentum 

Q 



PaUa + Pa U^a " 



dt 



PcUc + Pc 



dt 



Uxa ~l~ Ua 



df 



Pa + PaUl + 2paUa U,a 



dx 
dl 
dt 



dl 

dx ' 



Pc + PcUt + 2pJJc u 



One can find similar equations in Refs. 



20 



dl 
dt 



(38) 
(39) 
(40) 



26|. At this point one can check that Eqs. 



(!38|) - (HOl) are not sufficient to solve the problem. Indeed, we have only 3 equations for 4 
unknown values: the mode amplitudes Uza-, Uzs, Uzv, and the front perturbation /. Here we 
face the deficit of boundary conditions at the ablation front known in the problem of the RT 
instability in inertial confined fusion j?, 8]: in order to solve the problem we need one more 
condition at the front. Strictly speaking, one faces the same problem in the flame stability 
theory. Still, in the classical DL case of incompressible flow the missing condition may be 
easily guessed: this is the condition of a constant velocity of flame propagation with respect 
to the fuel mixture 

df 



u. 



Ref. 



dt 



0. 



(41) 



20| assumed a condition identical to (HTi) for a flame in a compressible flow. Surpris- 



ingly, that was not the only assumption of paper [20|. In addition to (14T!) . Ref. |20| assumed 
also a similar condition for the hot gas 

df 



u. 



dt 



0. 



(42) 
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In the limit of incompressible flow Eq. (l38l) reduces to 

Pa \Uza - j = Pc I Uzc ' 



dl 
dt 



(43) 



so that Eq. (1421) follows from fHT]) . However, in the case of compressible flow, these two 
equations are different. As a result, paper [20 1 suggested 5 equations for 4 unknowns, which 



made the problem over-defined. The paper 



26| criticized the solution [20], but it did not 



avoid an arbitrary assumption either. The paper [26| complemented the system ( !38|) - (HOj) 
by the perturbed equation of energy conservation, and the system of modes ( !3TI) - ( !36l) by 
the entropy mode (1271) . As a result, the approach of Ref. [261] involved 4 equations and 5 
unknowns. In order to obtain the solution, one mode was neglected in Ref. 26|]; namely, 
the entropy mode. That was just another arbitrary assumption; no physical explanation 
was provided in Ref. 26|] for this step. Thus, we come to the problem of correct additional 
boundary condition at a discontinuous deflagration front, which should replace fHTl) in the 
case of compressible flow. 

The missing equation may be derived only taking into account the perturbed equation 
of energy transfer, Eq. (|T71) . The term k'^KT in Eq. (|T71) may be neglected since it is small 
as kLc << 1 in the approximation of a discontinuous front. We split other terms into two 
groups: 



and 



~dS dS rf2 ^ ~ 

dz dz dz"^ 



apTS + {puz + pUz)T 



dp ^ 



dS 
dz 



dT ' 



We notice that perturbations in the form of front bending 



dz 



P 



-.dp 



MS 



^ dz^ ^ dz 



(44) 



(45) 



(46) 



turn the group Eq. fl44l) to zero. To find this, we take z-derivative of Eq. f[TT|) for the 
stationary deflagration front. We also notice that perturbations of density and temperature 
in the sound modes ( 13T1) - ( IMl) make the next order corrections in kL^ << 1 to the solution 
( l46l) . We substitute (H6l) into the perturbed continuity equation Eq. dHl), integrate 
over the deflagration front for zero-order terms in kL^ << 1 , and find 



{pU^ + pU^) - {pU^ + pU^)a = Crf{p - Pa), 



(47) 



11 



which corresponds to the perturbed equation of mass conservation at the front, Eq. (l38l) . 
In a similar way, substituting Eq. fj46|) into (|T5l) . f|T6l) and integrating, we may obtain the 
perturbed conservation laws fl39p . (1401) . As a next step we have to check, if the solution (1461) . 
fH7|) turn the second group, Eq. fHS]) . to zero. Substituting the solution fHBl) . fHTl) into fHS]) 
we find 

apTS + (pu^ + pUz)T^ = [{puz + puz)a - ofpa] (48) 
This combination turns to zero under the condition 

{pUz + pU^)a - (jfPa = 0. (49) 

Thus, solution in the form of front bending (146|) satisfies the equation of energy transfer 
with the additional condition 

PaUa + Pa(u,a-^^ =0. (50) 

Equation fl50|) plays the role of the DL condition modified for a compressible flow. In the 
incompressible limit. Ma << 1, density perturbations are negligible and Eq. (!50l) goes over 
to the classical DL condition Eq. ( HTl) . The physical meaning of Eq. ( 150|) is zero perturbation 
of the mass flux, which is mass burning rate in combustion or ablation rate in laser fusion. 
It is interesting that the obtained additional condition, Eq. (!50|) . does not depend on the 
type of energy release, see the derivation. We also stress that Eq. fl50l) is not a conservation 
law, and it should not be confused with energy conservation, though it was obtained from 
the equation of energy transfer. On the contrary, Eq. ( l50l) follows from the equation of 
energy transfer as an eigenvalue, similar to the incompressible case jsl, or even to the planar 
stationary Zeldovich - Frank-Kamenetski solution This result should be expected, since 
the deflagration speed (mass burning rate) is an eigenvalue of the basic equations ([1]) - ([3]), 
which does not follow from the conservation laws [1] . By this reason, the deflagration speed 
has to be obtained as an eigenvalue both for the unperturbed solution and for perturbations. 
This problem is intrinsic for deflagration, which is a discontinuity separating two subsonic 
flows. The problem of additional conditions does not arise in the case of detonation, for 
which one of the flows is supersonic. Unlike deflagration, dynamics of a detonation front is 
determined completely by the conservation laws. Thus, the reasoning of Ref. [26|] should 
work for detonation, but not for deflagration. The condition Eq. (1501) leads to an interesting 
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consequence. Substituting (150|) into the equation of dynamical pressure balance, Eq. (HOl) . 
we obtain 

P,{l-Mal)=P,{l~Mal). (51) 

Equation (ISTIl indicates the critical role of flow compression in the extreme CJ deflagration 
regime with Ma^ = 1. In the CJ regime, sound in the hot plasma propagates to the front 
as fast as it is drifted from the front by the flow. As a result, pressure perturbations in the 
cold plasma cannot be balanced by similar perturbations in the hot plasma. So, one should 
expect that the CJ regime is a critical one for the DL instability. For comparison, assuming 



extra condition fHTl) like in 20|, l25|], one obtains "+" instead of "— " in both sides of Eq. 



(15T|) . which means nothing special for the CJ deflagration. The special role of the CJ regime 



has been also indicated in 



26|. 



C. Solution to the stability problem 

The modes in the uniform flows (!32l) . (!33l) . (|37l) and the matching conditions at the 
deflagration front fl55]) . fHUl) . fISUl) determine the solution to the spectral problem for T = 
a/Uak. After heavy but straightforward algebra presented in the Appendix one can reduce 
this set of equations to 

r 0-r 



[Tr^s - e)(r + r^,) + (r - er^.) ^1 + - —^)\ = 0' (52) 

where the following designations are introduced 

Va = ^1 + Ma2(r2-1), (53) 



r^, = ^ll + Mal(^^-iy (54) 

In the limit of small Mach number Eq. ( l52ll reduces to the classical DL result. Numerical 
solution to (l52l) is shown in Fig. [2] versus the Mach number for different values of the expan- 
sion factor O = 4, 6,8 and the adiabatic exponent 7 = 5/3. For all expansion coefficients 
the instability growth rate decreases monotonically with the Mach number and it turns to 
zero in the adiabatic CJ regime of extremely fast deflagration at Ma^ = 1. Laser ablation 
in inertial confined fusion corresponds to the isothermal CJ regime with Ma^ = 1/7. In 
that case the DL instability is non-zero, but it is considerably weaker than in the classi- 
cal incompressible case. This is different from the previous results of Ref. [20t] predicting 
stronger DL instability at a discontinuous combustion front in a compressible flow. This 
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is also different from the results of Ref. 



2J] predicting no influence of compression effects 



for a discontinuous front of laser deflagration/ablation. According to Ref. 2J], compression 
effects come to play only as the second-order terms proportional to the finite deflagration 
thickness. On the contrary, we obtained decrease of the DL instability growth rate because 
of plasma compression already for a discontinuous deflagration front. Taking into account 
finite deflagration thickness we will, presumably, find strong influence of compression effects 
too. 

The numerical studies 



31 



32] of the D 



comment. The numerical data of Refs. 



31 



instability of a fast flame require a separate 



32[ demonstrated increase of the DL instability 



because of the compression effects for a fast flame of finite thickness with the Arrhenius 



reaction of the first order. Comparing 



31 



;he present theory to the numerical results of Refs. 



32l | one should notice that papers 



31 



32( 1 used the heating factor (temperature ratio) 



6fc = Tc/Ta as a fixed parameter instead of the expansion factor (density ratio) 6 = Pa/Pc 
used in the present paper. The heating factor is related to the expansion factor as 

76(9 - l)Mal 

e + 7(e-i)Mar ^ ^ 

The parameters G;, and G coincide for an incompressible flow, but they differ considerably 
when compression is strong. Choosing fixed G;, or G have quite different physical meaning, 
and we have to decide which parameter is more appropriate for the investigation. In a 
compressible flow, the deposited energy goes partly into the thermal energy (heating) and 
partly into the kinetic energy of the flow, which makes temperature ratio Gf, smaller than 
the density ratio G, see Eq. fl55|) . Still, the DL instability develops because of the density 
ratio: strong heating without expansion cannot produce the DL instability. Respectively, 
trying to keep the heating factor G5 fixed in the studies at large values of the Mach number, 
we inevitably have to increase the expansion factor G. Larger densities ratio G produces 
stronger DL instability, which makes a deceitful impression that compression effects enhance 
the DL instability. Such an approach may be, probably, justified for the combustion studies, 
but it definitely does not hold for laser ablation in inertial confined fusion. In the case of 
laser ablation, the density ratio G in plasma is specified by the frequency of laser light, 
which makes the condition of fixed expansion factor G quite natural. For comparison, we 
have plotted solution to the spectral problem, Eq. ( l52l) . for the fixed heating factor G;,. The 
respective plots are presented in Fig. [2] by the dashed lines. The dashed lines do show some 
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increase of the instability growth rate at small and moderate values of the Mach number 



similar to the numerical data of Refs. 



31 



32j. The numerical simulations Refs. 



31 



32] have 



been performed for combustion fronts of finite thickness. In order to compare the present 



case to Refs. 



31 



32| quantitatively, we also have to take into account finite thickness of the 



ablation region in inertial confined fusion. This is the subject for future work, which is in 
progress now. 

At the end of this subsection, we would like to point out that complete decrease of the 
DL instability in the CJ regime has an interesting physical interpretation. It is well-known 
that a deflagration front cannot propagate faster than the CJ deflagration On the other 
hand, it is also well-known that the DL instability increases the deflagration velocity at the 
nonlinear stage Then, what kind of nonlinear outcome we may expect for the DL 

instability at the CJ deflagration? The only possibility is that there is no DL instability in 
the CJ regime at all; we have obtained the same result solving the stability problem. 

D. Influence of source localization 

We have to make one more comment on the validity domain of the present results. In 
Sec. II we considered only localized energy sources Q^, which turn to zero in the hot plasma 
together with the first derivatives. This restriction is important for the present analysis. As 
a counter-example, let us consider a hypothetic energy source in the form 

nn = n,ip- p,y exp (^-Py^ ■ (56) 

Here pc is the final density of hot plasma, which is similar to density at the critical surface, 
the factor j3 » 1 plays the same role as Zeldovich number (scaled activation energy) in an 
Arrhenius reaction, u imitates the reaction order and Qa is some constant. The function (|56ll 
is constructed taking into account similarities with the Arrhenius reaction. In the case of the 
first-order reaction u = 1 and large Zeldovich number (3 » 1 the function (1561) is localized, 
but the first derivative of (156!) is not. Obviously, dQ^/dp is non-zero in the hot plasma 
with p = Pc for u = 1, which modifies the perturbed equation flM|) . As a consequence, 
the structure of sound modes in the hot plasma is modified too. In that case, within the 
discontinuity approach of kLc << 1, the leading term in ( l24l) is {OVLr/ dp)p = 0. This also 
leads to Pc = in the hot plasma instead of the adiabatic relation ( l26l) . We also obtain 
Pc = —k instead of fl5Ul) : and the structure of the acoustic wave in the hot plasma coincides 
with the incompressible case. The instability growth rate in that case is almost the same as 
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in the incompressible limit: the difference is only in the small terms Ma^ << 1. As we can 
see, taking the energy source in the form fl56|) with u = 1 we come to the results close to 



and different from the main conclusions of the present paper. Unfortunately, the 
work [2J| did not specify details of the energy release in the analysis. As another example, 
we can point out the numerical solution Refs. j3l|, 1321 . The solution of Ref . 3l|, |32| involved 



energy release described by an additional differential equation corresponding to a hypothetic 
first-order Arrhenius reaction. Similar to fl56l) with z/ = 1, perturbations of the energy source 



in Refs. 



31 



32l ] were not localized within the deflagration front. On the contrary, taking 



the second order u = 2 in Eq. ( l56i) instead of the first one, we have the energy release 
strongly localized together with its derivatives. For u > 1 in Eq. (1561) we recover the main 
results of the present analysis. Therefore, solution to the stability problem depends on the 
type of energy release. If the energy release is localized, then we obtain reduction of the 
DL instability by the compression effects shown in Fig. [2l If the source is not sufficiently 
localized, then we may obtain increase or decrease of the instability growth rate, or no effect 



at all, similar to the previous works 20|, l24, IMI, l32]- This reasoning is of special concern 
for the work [26], where the additional condition was not formulated as an equation at the 
front; but it took the form of an extra condition imposed in the flow of hot gas (plasnia). 
One cannot rule out, that the assumption of zero entropy perturbations adopted in Ref. [26 1 
holds for a certain type of energy source. 

A referee of the present work suggested also considering the case of isothermal sound in 
the uniform flow of hot plasma, Tc = 0, which may be obtained by appropriate modifications 
of the energy source. Such a condition is in line with the idea of isothermal plasma corona for 

2^. In that case Eq. (126|) should be replaced by P = {P/p)p 
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a planar ablation flow 

in the hot plasma and Ma^ by •yMal in Eqs. fl30l) . flSTl) and fl5^ . Respective solution to the 
stability problem is shown in Fig. [31 In that case complete stabilization of the DL instability 
happens already in the isothermal CJ regime, which is the hydrodynamic regime of laser 
ablation. We stress that such a result is not general, but it holds only for a specific type 
of energy source. Any attempt to obtain the same result in general by appealing to strong 
thermal conduction in plasma corona is incorrect as long as one works within the model 
of a discontinuous front. We remind that the condition of strong thermal conduction for 
perturbations implies kL^ » 1, which is just the opposite to the condition of a discontinuous 
front kLc << 1. Taking into account finite deflagration thickness and thermal conduction 
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like in ^, we can evaluate the cut-off wavelength of the ablation DL instability to be an 
order of magnitude larger than Lc- Therefore, the DL instability develops on length scales 
much larger than the distance from the ablation surface to the critical one. At such large 
length scales one may not treat sound perturbations as isothermal ones. 

Thus we come to the question if the energy source is sufficiently localized in laser fusion. 
The spacial damping rate of the laser light intensity is 33| 

with the absorption coefficient 

71^ n 



where n is electron concentration and ric is the critical value, for which plasma frequency is 
equal to the laser frequency 

Mark that laser radiation in ( 1571) propagates in the negative direction. Electron concen- 
tration in plasma is proportional to density n = p/M, where M is mass per one electron. 
Localization of energy source in (!57|) is due to the singularity of the damping rate (l58l) at 
n = ric, p = Pc- If we consider derivatives of the damping rate 

dK 3K 



dT 2T 



(60) 



^ - ^ + ^ (61) 

we find localization as strong as in (l58l) or even stronger. Thus, energy source in laser 
ablation is indeed strongly localized together with its derivatives, and the results of Fig. [2] 
hold for laser ablation. 



IV. SUMMARY 

In the present work we have solved the problem of the DL instability at a fast defla- 
gration in a compressible flow. The solution is obtained within the traditional model of a 
discontinuous front. Still, the approach of a discontinuous deflagration front suffers from 
the deficit of boundary conditions. We derive the missing condition as an eigenvalue of the 
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FIG. 1: a) Profiles of scaled density p/pa, tenl^rature T/Tc and energy release Q.lhQ,/^rnax for 
a defiagration front with = 10, k oc T^/^, and Ma^ = 1/7 (the isothermal CJ regime), b) The 




FIG. 2: Scaled instability growth rate F = a/Uak versus the Mach number in hot plasma for the 
fixed expansion factors 6 = 4, 6, 8. The dashed lines show similar solutions for the fixed heating 
factors @b = 4, 6, 8. The dashed-dotted line indicates the isothermal C J regime of laser ablation. 

equation of energy transfer. The derived condition corresponds to a constant mass flux at 
the deflagration front. In the limit of incompressible flow it goes over to the classical DL 
condition. We demonstrate that solution to the problem depends on the type of energy 
source. In the common case of a strongly localized source, compression effects make the DL 
instability considerably weaker. In particular, the DL instability growth rate is reduced for 
laser ablation in comparison with the classical incompressible case. The instability disap- 
pears completely in the Chapman- Jouguet regime of deflagration. However, if the energy 
source is not sufficiently localized, then it may influence properties of the sound perturba- 
tions in the hot plasma behind the front. In that case the properties of the instability also 
depend on a particular type of the energy source. 
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FIG. 3: Scaled instability growth rate F = a/Uak versus the Mach number in hot plasma for the 
fixed expansion factors = 4, 6, 8 in the case of adiabatic (dashed) and isothermal (solid) sound 
perturbations behind the front. The dashed-dotted line indicates the isothermal CJ regime of laser 
ablation. 
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APPENDIX: DERIVATION OF EQ. (f52]> . 

Here we present short derivation of Eq. (1521) . The spectral problem is determined by the 
equations (EOl) - (I32l), (JST]) - (iO]), ([501). Equation ^ may be also replaced by (EI]), which 
is more concise. Equations (1381) . (!50|) lead to 

UaPalc], + Pa{u,a " TUakf) = (Al) 
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f/,Pe/cL + Pc{u,c - rUakf) = (A2) 

We perform the following steps in calculations: 

1) Pc is expressed through Pa and Uza using ( 13T1) . ( 15T|) : 

2) fc/ is coupled to using (1A1I) . ( pTl) : 

3) using kf and ( 1A2I) . we relate Uzc and -Uza; 

4) using ([39]), ([32]), ([Ml), dSQ]), we relate u^c and m^^; 

5) everything is substituted into (137)) . 
Steps 1) - 4) in the calculations lead to 

Pc = aPa = -apPaUaUza, (A3) 

rUakf = (1 - Mal(3)uza, (A4) 

= (1 - Ma^/? + a/3Mal)uza, (A5) 

Mxc = tuza (— - ^^(1 - MalP)] , (A6) 
where the following designations have been introduced for a and (3: 

I -Mai ^ rfc ^ 
« = 1 = — + 1- A7 

Substituting (1A3P - (1A6P into Eq. (137)) we obtain solution to the spectral problem in the form 
of a single algebraic equation 

(A8) 

This equation may be further simplified. We introduce auxiliary designations: 
which leads to 

A : : lV= + Q = r^Mag- err/, + 9^(1 -Mag) 



Substituting ( 1A9I) -( 1ATI1) into ( 1A8I) and multiplying by (F — Or/c) we obtain 
af3{Vr]c - e)(l - Mag) + (F - Or^J [— + (1 - Mag/3) f ^ ® " ^ 



(A12) 
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Taking into account that 



k T + r]a 



we come to the final equation (152|) 



;i - Mal(3) 



(A13) 



l + Va 



r _ e - 1 
e F~ 



(AM) 



This gives solution to the spectral problem in a form of one concise algebraic equation. 
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